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Abstract 

A  state  of  charge  (SOC)  algorithm  is  described  and  implemented  for  a  nickel  metal  hydride  (NiMH)  battery  system.  The  essential  elements 
of  the  algorithm  were  first  used  as  part  of  the  GM  Precept  hybrid  electric  vehicle  program,  which  culminated  in  an  80-mile-per-gallon, 
5-passenger,  technology-demonstration  vehicle.  The  algorithm  is  based  on  a  simple  equivalent-circuit  representation  of  the  battery  pack, 
with  the  parameters  adapted  by  means  of  on-line  least-square  regression.  The  output  of  the  algorithm  is  the  SOC,  which  is  a  measure  of 
the  available  energy  within  the  battery  system,  and  the  power  capability  for  subsequent  charge  or  discharge,  which  enables  efficient  hybrid 
vehicle  operation  through  optimal  torque  allocation.  The  behavior  of  the  algorithm  in  terms  of  convergence,  accuracy,  and  robustness  is 
examined;  for  these  analyses,  the  power  required  from  the  battery  for  the  Precept  to  complete  a  variety  of  drive  schedules  provides  the 
power  versus  time  trace.  Voltage  hysteresis  is  shown  to  play  a  critical  role,  and  this  effect  makes  determining  the  SOC  particularly  difficult 
for  NiMH  batteries  relative  to  other  battery  systems  (e.g.,  lead  acid  and  lithium  ion). 
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1.  Introduction 

There  are  two  terms  of  immediate  interest  to  the  user 
community  of  propulsion  battery  systems:  state  of  charge 
(SOC)  and  state  of  health  (SOH).  The  SOC  corresponds  to 
the  stored  charge  available  to  do  work  relative  to  that  which 
is  available  after  the  battery  has  been  fully  charged;  this  def¬ 
inition  is  made  precise  in  the  model  formulation  to  follow. 
SOC  can  be  viewed  as  a  thermodynamic  quantity,  enabling 
one  to  assess  the  potential  energy  of  the  system.  SOH  is 
a  term  that  is  becoming  more  commonly  used  within  the 
battery  community,  but  which  has  not  to  date  been  clearly 
defined.  Generally,  SOH  is  used  to  imply  that  one  can  de¬ 
duce  how  well  the  battery  system  is  functioning  relative  to 
its  nominal  (rated)  and  end  (failed)  states.  For  our  purposes, 
we  assume  that  we  can  represent  the  SOH  if  we  have  a 
method  to  identify  the  impedance  spectrum  for  the  battery 
system  over  the  frequency  range  of  interest  in  an  on-line 
(adaptive)  manner.  Hence,  knowing  the  change  in  the  SOH 
with  time  may  be  viewed  as  enabling  one  to  assess  the  in¬ 
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crease  in  irreversible  losses  that  is  inherent  in  the  aging  of 
batteries.  SOH  is  relevant  to  on-board  diagnostics.  Thus,  the 
system  thermodynamics  allow  one  to  assess  the  potential 
energy  of  the  system,  and  irreversible  losses  can  be  assessed 
once  the  impedance  spectrum  is  known.  The  parameters  that 
are  used  to  characterize  the  impedance  spectrum  can  be  re¬ 
gressed  by  means  of  a  system  identification  scheme.  The 
algorithm  developed  in  this  work  was  applied  to  the  Precept 
hybrid-electric- vehicle  battery  system.  (See  Figs.  1  and  2 
for  the  vehicle  and  NiMH  battery  pack,  respectively,  [1].) 
The  overall  approach  should  also  be  useful  for  lead  acid  and 
lithium  ion  systems  as  well.  For  an  analysis,  overview,  and 
summary  of  challenges  associated  with  battery  monitoring 
and  energy  management  for  electrical  power  systems,  the 
paper  by  Meissner  and  Richter  is  quite  helpful  [2]. 

We  shall  determine  the  battery’s  SOC  and  SOH  (the  lat¬ 
ter  insofar  as  the  high-frequency  resistance  serves  a  proxy 
for  the  SOH)  based  on  least-squares  recursive  identification 
techniques.  The  equivalent  circuit  shown  in  Fig.  4  shall  be 
employed  to  represent  the  essential  features  of  the  battery 
system.  Reviews  on  the  substantially  interdependent  fields 
of  recursive  identification,  adaptive  filters,  optimal  estima¬ 
tors,  and  model-reference  adaptive  systems  (MRAS’s)  can 
be  found  in  references  [3-9].  An  objective  of  this  work 
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Fig.  1.  Rendering  of  the  Precept  hybrid  electric.  The  battery  system  is  located  under  the  front  seats.  The  front  axel  has  a  coaxial  electric  traction  system, 
and  the  rear  axel  receives  torque  through  a  traction  system  comprising  an  electric  motor,  a  compression  ignition  heat  engine,  and  a  manual  transmission 
that  is  automatically  shifted  (actuated).  Regenerative  braking  takes  place  at  both  axels. 


was  to  develop  an  equivalent-circuit  model  that  can  be  used 
with  any  recursive  identification  scheme;  the  model  is  then 
combined  with  the  least- squares  method  so  as  to  regress 
the  relevant  parameters  consistent  with  a  model  reference 
adaptive  system  and  application  in  an  embedded  controller 
[7].  More  complex  models  that  include  the  numerical  so¬ 
lution  to  coupled  partial  differential  equations  for  the  pur¬ 


poses  of  describing  thermodynamic,  kinetic,  and  transport 
phenomena  from  a  microscopic  basis  [10,11]  would  likely 
lead  to  a  prohibitively  complex  MRAS,  as  affordable  em¬ 
bedded  controllers  do  not  have  the  necessary  storage  and 
execution-time  capabilities  today  (or  in  the  immediate  fu¬ 
ture),  although  promising  work  in  this  area  is  taking  place 
[12].  The  least-squares  formulation  of  this  study  amounts  to 


Fig.  2.  Battery  pack  system.  The  28  battery  modules  are  liquid  cooled;  the  coolant  enters  and  exits  at  the  rear  of  the  pack.  Controllers,  high-voltage 
devices  (contactors,  relays,  etc.)  are  located  in  front  and  behind  the  modules. 
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Percent  State  of  Discharge 


Fig.  3.  Open-circuit  voltages.  The  upper  panel  illustrates  the  flatness  of  the  open-circuit  voltage  vs.  SOC.  The  lower  plot  utilizes  an  expanded  ordinate. 


a  fixed  frame  of  data  points,  and  only  two  parameters  are  re¬ 
gressed  (the  high-frequency  resistance  and  the  open-circuit 
voltage).  As  will  be  discussed  in  detail,  the  NiMH  sys¬ 
tem  is  complicated  by  hysteresis  in  the  open-circuit  volt¬ 
age  (cf.  Fig.  3),  [13-17]  which  must  also  be  regressed 
adaptively. 

2.  Model  formulation 

The  model  used  for  this  estimation  system  consists 
of  two  parts.  There  is  an  electrical  circuit  model  that  is 
used  to  describe  the  relationship  between  the  currents  and 
voltages  observed  at  the  terminals  of  the  battery  and  a 
coulomb-accumulation  model  that  describes  the  open-circuit 
voltage  based  on  the  history  of  currents  seen  by  the  battery, 
including  self  discharge  and  current  inefficiency  on  charge. 
The  electrical  circuit  model  is  illustrated  in  Fig.  4.  Because 


the  SOC  is  a  function  of  the  open-circuit  potential,  extrac¬ 
tion  of  the  open-circuit  voltage  by  least-squares  regression 
leads  to  an  estimate  of  the  SOC.  The  two  SOC  values 
are  referred  to  as  SOCc  (for  the  coulomb-accumulation 


Fig.  4.  Equivalent  circuit  for  the  battery  system.  The  representation 
does  not  include  self  discharge,  current  inefficiency,  temperature  and 
SOC  dependence  of  “circuit  elements”,  and  hysteresis  in  open-circuit 
potential. 
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contribution)  and  SOCy  (for  the  voltage-based  contribu¬ 
tion).  Both  contributions  yield  useful  information  regarding 
the  SOC;  thus  a  composite  SOC  is  calculated.  To  find  this 
composite  SOC,  the  following  is  utilized: 

SOC  =  w( SOCc)  +  (1  -  w)( SOCy),  (1) 

where  w  is  a  weighting  factor.  For  SOCc  [18], 

SOCc(0  =  SOC(f-  AO 


Jt-At 

Time  is  represented  by  t,  and  I  denotes  current;  discharge 
currents  are  taken  as  positive.  It  is  important  to  note  that 
the  first  term  on  the  right  is  SOC,  not  SOCc;  the  use  of 
the  composite  SOC  as  the  first  term  on  the  right  strongly 
couples  the  voltage-based  and  coulomb-accumulation  mod¬ 
els.  For  the  module  of  this  example,  the  nominal  capac¬ 
ity  Ahnominai  is  12.5  Ah  (35  °C,  C/3  discharge  rate;  this 
quantity  will  have  a  mild  temperature  dependence  over  the 
range  of  application),  and  the  self-discharge  rate  is  described 
by  [18] 

SD  =  koexp  (-yy)  SOC, 


100  .. ??//  +  sd(T) 

Annominai 


d  t 

3600’ 


(2) 


Rct  [19].  An  analysis  of  proton  transport  in  nickel  hydroxide 
can  be  found  in  the  treatment  by  Weidner  and  Timmermann 
[24]. 

It  is  necessary  to  have  a  default  value  for  the  resistance; 
the  approach  for  extracting  the  resistance  from  recorded  data 
will  be  discussed  below.  For  the  NiMH  battery  example 
application  in  this  work,  the  resistance  data  were  fit  to  the 
following  equation  [18]: 


^module  —  ^  ^  Clj  SOC^ 

The  pack  resistance  is  given  by  the  number  of  modules 
multiplied  by  /^module  as  other  resistances  (e.g.,  harnesses, 
etc.)  have  been  allocated  to  the  module  resistance.  (For  the 
pack  utilized  in  this  example,  28  modules  were  employed. 
Thus,  R  =  A^mociuies ^modules ?  where  Amoc iuies  =  28.)  The 
activation  energy  Ea,r  corresponds  to  —6,017  cal  mol-1.  It 
is  convenient  to  choose  the  gas  constant  Rg  with  consistent 
units,  1.987  cal  mol-1  K-1).  The  reference  temperature  Tre f 
is  chosen  as  35  °C,  or  308  K.  The  coefficients  in  the  sum¬ 
mation  are  given  below.  For  all  of  the  results  shown  in  this 
work,  the  regressed  resistance  was  subsequently  bounded 
Such  that  0.5A7?module  5:  R  —  lAV7?module* 


exp 


Ea,r  (  1 
'  Re  \T 


1 

^ref 


with  ko  =  1.0683  x  107  per  hour  and  EA,s/Rg  =  6,  789  K. 
The  current  efficiency  77/  is  effectively  unity  for  the  con¬ 
ditions  near  50%  SOC,  as  confirmed  by  measurements;  a 
value  of  1  is  used  for  discharge,  and  0.99  for  charge;  current 
efficiencies  are  expected  to  have  a  mild  temperature  depen¬ 
dence  (for  the  conditions  of  this  work)  and  decrease  with 
increasing  SOC’s.  The  factor  3600  has  units  of  s/h,  and  the 
factor  100  is  employed  to  keep  a  consistent  percent  basis. 

To  extract  the  voltage-based  SOCy,  the  solution  for  the 
voltage  associated  with  the  equivalent  circuit  shown  in  Fig.  4 
is  employed  [18-20]: 


V  =  Vo— IR  +  A 


C=o 


/ 

h=t 


/(?)exp[-5(f-0]d£ 


(3) 


The  first  two  terms  on  the  right  side  give  rise  to  an  ohmic 
description  of  the  battery,  as  the  battery  voltage  V  is  related 
to  the  open-circuit  potential  Vo  reduced  by  the  ohmic  drop 
IR,  where  R  is  the  battery  resistance.  The  last  term  on  the 
right  side  corresponds  to  a  superposition  integral,  through 
which  past  currents  influence  the  open-circuit  potential  be¬ 
yond  the  first-order  effect  of  changing  the  average  state  of 
charge  characterizing  the  electrodes.  (Because  of  the  expo¬ 
nential  weighting  function,  the  influence  of  early  currents  is 
less  than  that  of  more  recent  currents.)  An  incremental  im¬ 
provement  in  the  generality  of  the  equivalent  circuit  would 
be  to  add  a  capacitor  in  series  with  the  interfacial  resistor  Rct 
or  a  frequency-dependent  Warburg  impedance.  Such  addi¬ 
tions  are  intended  to  approximate  mass-transfer  impedance 
[21-23];  the  basic  form  of  Eq.  (3)  does  not  change  substan¬ 
tially  when  a  capacitor  is  added  in  series  with  the  resistor 


j  cij  (ohms) 


0 

4.1252  x  10“2 

1 

8.9691  x  1(T4 

2 

1.6760  x  10“5 

3 

-1.4435  x  10“7 

4 

4.7223  x  IQ-10 

The  open-circuit  voltage  (OCV)  is  described  by  a  modi¬ 
fied  Nernst  equation,  with  the  inclusion  of  the  influence  of 
entropy  of  reaction  and  an  empirical  expression  to  capture 
the  salient  features  associated  with  voltage  hysteresis: 


Voc,cell  =  Function (T,  SOC,  VH) 


=  U°+VH  + 


RgTreiln  SOC  -  77 
neF  ln  100 -SOC 


-  (7)  (SOC) 


+  (T-  Tief) 


AS 


nQF 


(4) 


Faraday’s  constant  (96,487  C  mol-1)  corresponds  to  F,  ne 
denotes  the  electrons  according  to  the  cell  reaction  = 
1  for  this  study),  and  AS  refers  to  the  entropy  of  the  cell 
reaction,  22  J  mol-1  K-1.  (For  consistent  units,  a  value  of 
8.314  J  mol-1  K-1  should  be  used  for  Rg  in  Eq.  (4).)  To 
construct  the  OCV  curve  shown  in  Fig.  3,  U°  =  1.37  V,  77  = 
8  (taken  as  dimensionless,  as  is  the  percent  state  of  charge), 
and  y  =  0.04  V.  The  number  of  cells  in  the  battery  pack  for 
this  example  is  given  by  N  =  280;  thus  the  pack  open-circuit 
voltage  is  given  by  V0c  =  Woc,cell-  The  last  term  in  Eq.  (4) 
corresponds  to  a  linearization  of  the  open-circuit  potential 
about  the  reference  temperature,  35  °C  (308  K)  in  this  work. 
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For  the  hysteresis  contribution,  we  construct  the  following 
first-order  differential  equation  to  formulate  a  cell  hysteresis 
voltage  Vr: 

^  =  -Piriil  -  sSD)  [yH,max  +  sign(/)Vk]  •  (5) 

This  equation  constructs  a  varying  hysteresis  voltage; 
for  this  work,  the  hysteresis  voltage  is  set  up  so  that  for 
prolonged  charge  currents,  or  short  but  very  large  charge 
currents,  the  hysteresis  voltage  tends  to  50  mV  per  cell 
(for  280  cells  per  pack,  this  corresponds  to  11.2  V);  thus 
^H,max  =  50  mV  per  cell.  The  exact  opposite  holds  for 
discharge  (positive)  currents.  Note  also  that  if  the  current 
remains  at  zero  for  a  long  time,  the  hysteresis  voltage 
tends  to  the  charge-decreasing  condition  (—50  mV  per  cell) 
through  self-discharge.  The  ±50  mV  per  cell  corresponds 
to  the  100  mV  per  cell  difference  between  the  upper  and 
lower  open-circuit  voltage  curves  displayed  in  Fig.  3;  the 
middle  curve  of  Fig.  3  corresponds  to  Vr  =  0. 

The  description  of  the  governing  equations  is  now  com¬ 
plete.  The  next  section  provides  details  regarding  implemen¬ 
tation  as  an  SOC  algorithm. 


3.  Implementation 


First,  the  coulomb  integration  expression  2  is  recast  as: 


SOCc(0  =  SOC (t  -  At) 

_  '10()  [(»?//, -i +>?//, )/2] 

Ahnomjna] 


±  ( Sv)t-At 


At 

3600’ 


The  difference  between  the  present  time  and  the  last 
recorded  time  is  given  by  At.  Next,  the  measured  voltage 
expression  is  recast  for  on-line  regression.  First,  it  is  help¬ 
ful  to  note  that  Eq.  (3)  can  be  rewritten  as  the  following 
recursion  relation  for  evaluation  purposes: 


V\t  =  (Voc-IR)t-(I\tA)At 

±  exp(—BAt)[V  -  ( Voc  -  IR)]t-At,  (6) 

where  the  subscripts  t  and  t  —  At  denote  the  time  at  which  the 
quantities  are  to  be  evaluated.  This  equation  is  a  particularly 
simple  recursion  relation  in  which  only  variables  calculated 
at  the  previous  time  step  are  required  to  calculate  the  voltage 
at  time  t.  The  equation  is  derived  in  the  Appendix  A.  To 
implement  Eq.  (6),  one  replaces  the  battery  voltages  with 
measured  values: 

yto  be  used  in  regression 

=  V  =  Vm™smed\t  +  (I\,A)At 

-  exp ( —  fi A 0 [ Vmeasured  -  (Voc  -  IR)\t_Al 
=  (V0c  -  iR)t 

Thus,  the  regression  analysis  to  determine  the  open-circuit 
potential  and  resistance  is  based  on  the  voltage  V,  the  re¬ 


gression  voltage,  and  a  least-squares  analysis  of  the  cor¬ 
rected  voltage  data  should  yield  a  good  approximation  for  the 
ohmic  resistance  and  open-circuit  potential.  For  all  results 
shown  in  this  work,  A  =  0.0229  F-1and  B  =  0.0366  s-1 
[18]. 

The  open-circuit  voltage  Voc  and  the  resistance  R  are 
found  using  a  least- squares  approach,  applied  to  data  corre¬ 
sponding  to  a  specified  time  interval.  The  following  defini¬ 
tions  are  applied  [25]: 


n 


Zb 


7=1 


1 

SJV  =  ~ 
n 


where  n  represents  the  number  of  recorded  current-potential 
data  points  to  be  included  in  the  extraction  of  the  open-circuit 
voltage  Voc  and  the  resistance  R.  For  all  analyses  presented 
in  this  work,  the  time  step  is  1  s  and  n  =  91.  The  summations 
are  made  effectively  recursive  by  recognizing  that  a  first-in, 
first-out  method  keeps  the  memory  full  of  the  most  recent  n 
data  points  and  does  not  require  a  full  summation  over  all 
components  within  the  sum  at  each  time  step.  Using  these 
expressions,  one  obtains  the  following, 


siv  -  SjSy 
SII  -  (si)2 


(7) 


and 


_  susy  -  .S'lV-n 

s\\  -  Vi)2 


(8) 


These  equations  fail  when  the  variance  sn  —  si  =  0,  or 
when  this  quantity  is  nearly  zero.  In  addition,  the  equations 
can  fail  to  provide  a  reasonable  result  when  the  many  of  the 
recorded  currents  used  in  the  least-squares  analysis  are  of 
similar  value  with  the  exception  of  one  or  two  data  points; 
thus  in  application  it  is  important  not  to  update  (regress) 
R  and  Voc  when  the  variance  is  small  or  there  is  a  large 
skewness  [26]  in  the  data. 

The  hysteresis  contribution  is  now  addressed.  The  follow¬ 
ing  formulation  is  applied: 


Vn(t)  =  u;H{VH-j8(AO(i7//-SD)[VHfniax+sign(7)VH]h-Ar 


+  (1  -  wh) 


Voc  -  U 


ref 


RgTref  SOC  -  77 
neF  ln  100 -SOC 


+  (y)S0C-(r-rref) 


AS 

~nF 


t—At 


(9) 
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The  subscript  on  the  large  brackets  and  braces  indicates 
that  values  to  the  right  of  this  equation  can  be  evaluated  at  the 
previous  time  step.  This  equation  is  not  a  straight-forward 
(numerical)  time  integration  of  Eq.  (5)  unless  the  weight¬ 
ing  factor  von  is  set  to  unity.  Thus  the  second  line  of  the 
above  equation  allows  for  a  correction  to  the  extraction  of 
the  hysteresis  voltage  through  the  recognition  that  the  pre¬ 
vious  time  step  value  for  the  SOC  can  be  used  to  calculate 
an  open-circuit  voltage  (cf.  Eq.  (4)) — this  back-calculated 
open-circuit  voltage  provides  an  adaptive  routine  for  deduc¬ 
ing  VH. 

The  next  step  is  to  transform  the  open-circuit  voltage  into 
a  voltage-based  state  of  charge,  including  the  hysteresis  volt¬ 
age.  To  do  this,  the  modified  Nernst  Eq.  (4)  is  inverted  to 
yield  the  voltage-based  SOC.  Rather  than  numerically  solve 
a  non-linear  equation  in  an  embedded  controller,  one  can 
simply  use  a  look-up  table  to  represent  the  non-hysteresis 
portion  of  the  modified  Nernst  equation. 

Both  the  combined  SOC  (Eq.  (1))  and  the  hysteresis  volt¬ 
age  (Eq.  (9))  utilize  weighting  factors  (w  and  wn,  respec¬ 
tively),  and  the  influence  of  the  time  step  size  needs  to 
be  addressed.  That  is,  if  very  small  time  steps  are  em¬ 
ployed,  then  the  weighting  factors  should  be  altered  so  that 
time-dependent  quantities  are  not  lost  from  the  calculation 
and  instead  are  allowed  to  evolve  in  accordance  with  their 
particular  time  constants.  The  following  approach  is  applied 
in  this  work: 


Fig.  5.  Voltage-based  portion  of  the  SOC  algorithm.  The  recursive  nature 
of  the  integration  is  such  that  only  quantities  time  steps  k  and  k—  1  are 
required.  The  error  is  minimized  by  least-squares  regression  of  Vqc  and  R. 


low  (e.g.,  less  than  0.1  A).  In  this  case,  the  voltage-based 
state  of  charge  is  likely  to  provide  the  most  accurate  es¬ 
timate  for  the  combined  SOC — an  extreme  case  would  be 
for  a  battery  sitting  at  rest  for  many  days,  in  which  case 
the  chemistry  is  dictated  by  self-discharge  processes,  and 
the  charge-decreasing  Nernst  equation  provides  an  accurate 
estimation  for  the  battery’s  open-circuit  potential.  In  such 
cases,  w  can  be  set  to  a  value  near  0.5,  which  is  what  we 
have  employed  as  a  lower  bound  on  w. 

The  quantity  aw  is  obtained  as  follows.  First,  we  note  that 

^Vqc  ^^y^g^ref 

3(SOC)  ~  ^  nQF 

Thus  we  construct  aw  as  follows: 


100 

SOC(100  -  SOC) 


W  =  Wmax  -  aw(At)  and  wH  =  WH,max  -  un(At). 

The  weighting  factors  are  bounded  between  0  and  1  here. 
The  factor  an  is  taken  as  0.005  s-1  for  all  of  the  plots  in  this 
work.  For  short  times  (i.e.,  before  the  regression  analysis  al¬ 
lows  for  accurate  fitting  of  the  resistance  and  open-circuit 
potential,  91  s  in  this  work),  w  is  set  to  unity,  and  the  SOC 
is  calculated  for  this  short  time  based  solely  on  coulomb  in¬ 
tegration  from  the  previously  stored  SOC.  Another  special 
case  results  when  the  battery  has  been  at  rest  for  prolonged 
periods.  In  this  case,  the  hysteresis  model  returns  a  value  of 
Vh  that  is  less  than  —45  mV  due  to  self  discharge  (when  s  is 
greater  than  zero),  and  the  magnitude  of  the  current  is  quite 


For  SOC  >  50%  :  a 


W  I50%SOC 


For  SOC  <  50%  :  min 


10, 


li 


4 1 50%: 


SOC  J 


xc*%l50%SOC 


The  quantity  10  in  this  equation  corresponds  to  an  upper 
bound  used  in  this  work,  and  aw  at  50%  SOC  was  set  to 
0.001  s-1.  This  approach  yields  larger  values  of  aw  as  the 
SOC  declines  below  50%,  consistent  with  placing  a  greater 
emphasis  on  the  voltage-base  SOC  when  the  OCV  is  more 
sensitive  to  changes  in  SOC. 

A  high  level  summary  of  the  algorithm  is  depicted  in 
Fig.  5  (for  the  more  complicated  voltage-based  SOC  asso¬ 
ciated  with  Eq.  (6))  and  Fig.  6  (for  the  entire  algorithm). 


Inputs  at 
time  tfr 

Ik,Vk,Tk,SOCk_x 

A 


socv 

1.  Extract  R  and  Voc  [requires 

parameter  that  are  functions 

of  (tt,Tk,SOCk_x)} 

2.  Determine  SOCv  from 

Voc=VH{Ik)  +  V0{SOCr) 

SOCc  by  coulomb  integration 
*  (including  self  discharge  and 
current  inefficiency) 


w(t,SOCkA,IbVk) 


■> 


\-w(t,SOCk_xJbVk) 


unit  delay  z  1 


> 


SOCk  (and  power 
capability) 


Fig.  6.  Schematic  overview  of  algorithm.  The  outputs  are  the  combined  SOC,  the  charge  and  discharge  power  capabilities,  w  is  a  function  of  t,  SOQt-i, 

4,  V*. 
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Similar  schematics  can  be  found  in  the  controls  literature 
for  algorithm  description  [3-9] . 


4.  Results  and  analysis 

For  all  of  the  plots  shown,  the  parameters  values  listed  in 
Table  1  were  employed  unless  otherwise  stated.  Before  ana¬ 
lyzing  the  specific  application  of  the  algorithm,  an  overview 
of  the  data  shown  in  Fig.  7  provides  a  compelling  reason  as 
to  why  this  approach  is  promising.  A  drive  schedule  similar 
to  that  used  to  determine  the  fuel  economy  values  for  the 
Precept  HEV  yielded  the  data  points  depicted  in  Fig.  7. 
Note  that  the  points  are  substantially  co-linear,  and  a  simple 
ohmic-battery  description  ( V  =  Voc  — IR )  fit  to  the  data  cap¬ 
tures  the  basic  trend  in  the  data.  This  result  that  is  partially 
enabled  by  the  small  variation  in  the  SOC  throughout  the 
test;  if  large  variations  in  SOC  are  manifest,  the  Voc  would 
not  be  substantially  constant,  and  the  ohmic-battery  descrip¬ 
tion  would  be  less  insightful.  Fig.  8  can  be  viewed  as  the  base 
case  for  this  work,  as  it  shows  the  results  for  the  drive  sched¬ 
ule  of  Fig.  7.  The  battery  was  first  charged  to  50%  SOC,  then 
the  battery  was  stimulated  using  a  power  signal  that  approxi¬ 
mated  the  operation  of  the  battery  when  the  vehicle  is  driven 
on  the  FTP  driving  schedule.  The  measured  current  and  volt¬ 
age  are  shown  in  the  upper  plot,  along  with  the  SOC  deduced 
by  purely  coulomb  counting  and  the  open  circuit  voltage 
extracted  from  the  data  using  the  full  algorithm.  It  should  be 


Table  1 

Quantities,  values,  and  units  unless  otherwise  specified 


Value 

Units 

Quantity 

50 

mV 

Half  the  maximum  cell  hysteresis,  V^max 

12.5 

Ah 

Nominal  capacity  at  Tref 

0.0229 

C_1 

A,  inverse  capacitance 

0.0366 

s-1 

B,  inverse  time  constant 

10 

Cells  per  module  ( N  =  280  cells  per  pack) 

28 

Modules  per  pack 

0.99 

Current  efficiency  at  Tref 

35 

°C 

Reference  temperature  Tref 

2.47  x  10-5 

C-1 

P  for  charge-augmented  operation 

3.70  x  10-5 

C-1 

p  for  charge-depleted  operation 

0 

Ah 

s,  self  discharge  multiplier  for  hysteresis 
expression 

0.650 

Q 

Default  high-frequency  pack  resistance  R 

0.001 

s-1 

aw  at  50%  SOC 

0.005 

s_1 

m 

stressed  that  we  do  not  have  a  precise  method  to  determine 
the  error  in  an  SOC  routine,  as  one  cannot  fully  discharge 
the  batteries  at  low  rates  without  damaging  them.  Our  expe¬ 
rience  is  that  for  the  conditions  of  this  work,  the  relatively 
short  times  (less  than  a  few  hours)  of  power  cycling  allows 
us  to  use  the  pure  coulomb  counting  as  a  sufficiently  accurate 
measure  of  the  SOC.  For  longer  durations  (e.g.,  days),  this 
would  not  be  the  case.  The  lower  plot  in  Fig.  8  provides  the 
error  based  on  the  assumption  that  the  pure  coulomb  count¬ 
ing  returns  the  correct  answer.  The  percent  error  is  taken  as 


Fig.  7.  Battery  system  current-voltage  relation  over  a  representative  drive  cycle.  The  pack  power  was  allocated  based  on  vehicle  modeling  for  the  Precept 
vehicle  over  a  drive  cycle  of  the  following  characteristics:  1384  s  city  schedule  (FUDS),  10  min  rest,  505  s  of  city  schedule  (FUDS),  and  two  highway 
schedules  (HWYFET).  The  collected  current  and  voltage  data  for  the  described  schedules  correspond  to  the  plotted  data  points.  The  fit  line  indicates  an 
open-circuit  voltage  of  387  V,  reflecting  about  50%  SOC. 
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Fig.  8.  Algorithm  applied  to  data  shown  in  Fig.  7.  The  voltage  and  current  correspond  to  pack  data;  with  the  exception  of  the  percent  error,  all  other 
results  are  outputs  of  the  algorithm. 


100(SOCpure  coulomb  counting  SOC)/SOCpure  oulomb  counting- 
Other  quantities  plotted  in  the  lower  portion  of  Fig.  8  include 
the  extracted  module  resistance  (about  23  m£2  per  module 
throughout  this  work,  corresponding  to  a  pack  resistance  of 
650  m£2  as  shown  in  Fig.  7),  the  regressed  combined  state 


of  charge  (SOC),  and  the  regressed  voltage-based  state  of 
charge  (SOCy).  The  hysteresis  voltage  declines  from  nearly 
50  mV  per  cell  initially,  reflecting  a  battery  that  had  just 
been  charged,  indicating  that  the  operation  of  the  battery 
system  is  largely  charge  sustaining  throughout  the  test. 
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Fig.  9.  Importance  of  hysteresis  voltage.  The  hysteresis  voltage  is  set  to  zero  throughout  the  algorithm  application;  all  other  conditions  are  identical  to 
Fig.  8.  The  percent  SOC  error  is  seen  to  be  much  larger  that  that  of  Fig.  8. 
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Fig.  10.  Effect  of  initial  hysteresis  voltage  value.  Although  the  battery  was  charged  immediately  prior  to  the  experiment,  the  hysteresis  voltage  initially 
was  set  to  zero  (cf.  Fig.  8). 
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The  importance  of  the  hysteresis  in  extracting  the  correct 
SOC  is  made  clear  in  Fig.  9,  where  Vh  has  been  set  to  zero. 
The  SOC  initially  is  correct,  as  it  was  set  to  the  correct  50% 
value,  but  by  the  end  of  the  run  it  is  seen  that  error  in  the  SOC 
is  nearly  50%.  The  initial  value  for  the  hysteresis  voltage  is 
not  critical,  and  the  adaptability  of  the  algorithm  in  terms  of 


extracting  Vh  is  demonstrated  in  Fig.  10.  We  see  that  when 
the  initial  value  for  Vh  is  set  to  0  mV  per  cell,  instead  of  the 
correct  value  of  +50  mV  per  cell  after  charge,  the  algorithm 
tends  to  “converge”  in  a  couple  minutes,  and  the  error  in 
the  SOC  is  never  prohibitively  large.  This  implies  a  rather 
stable  algorithm. 


Time,  s 


Fig.  11.  Robustness  test.  As  in  the  previous  figure,  the  hysteresis  voltage  was  set  to  zero  initially.  In  addition,  the  beginning  SOC  used  to  start  the 
algorithm  was  set  to  10%  (instead  of  the  actual  50%).  Upper  figure:  aw  =  0.001s_1.  Lower  figure:  aw  =  0.01s-1  (more  adaptive).  The  error  in  the 
regressed  SOC  is  less  than  10%  within  100  s. 
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A  more  severe  test  of  the  algorithm  rate  of  convergence 
and  stability  is  depicted  in  Fig.  11,  where  the  initial  SOC  for 
the  algorithm  is  set  to  10%,  even  though  the  actual  initial 
SOC  is  50%.  By  comparing  the  upper  and  lower  plots,  we 
see  that  in  order  to  speed  convergence  of  the  algorithm  to  less 
than  a  couple  minutes,  the  weighting  of  the  voltage-based 


portion  of  the  algorithm  must  be  increased.  The  lower  panel 
demonstrates  superior  convergence  properties  with  the  larger 
value  of  aWj 5o%soc  =  0.01s-1.  However,  the  combined 
SOC  shows  more  variation  (albeit  still  acceptably  stable), 
and  an  expected  tradeoff  between  stability  and  adaptability 
must  be  comprehended  in  application.  We  have  found  it 
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Fig.  12.  Effect  of  prior  discharge  on  algorithm.  After  the  nearly  2000  s  discharge,  the  power  profile  is  very  similar  to  that  of  Fig.  8. 
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useful  to  leave  au;,50%soc  at  0.001  s-1  unless  a  need  arises 
to  reset  the  algorithm  (e.g.,  a  controller  reset  command  as 
might  be  expected  in  a  service  environment). 

The  effect  of  prior  discharge  on  the  algorithm  behavior, 
unlike  the  early  charge  event  investigated  in  Fig.  8,  is  ex¬ 
plored  in  Fig.  12.  After  the  nearly  2000  s  discharge,  the 


power  profile  of  Fig.  12  is  very  similar  to  that  of  Fig.  8.  In 
general,  the  results  are  similar  to  those  of  Fig.  8,  although 
the  error  is  a  bit  larger.  Further  tuning  of  the  algorithm  [4] 
can  reduce  this  error  (e.g.,  deceasing  ^discharge  reduces  the 
error  in  this  case).  The  last  investigation  involves  applica¬ 
tion  of  the  algorithm  to  a  high  power  profile,  characteristic 
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Fig.  13.  Algorithm  applied  to  high  power  schedule.  The  power  after  charge  (i.e.,  after  3500  s)  corresponds  to  a  US06  schedule  for  the  Precept. 
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of  what  the  Precept  sees  over  the  aggressive  US06  sched¬ 
ule.  The  magnitude  of  the  current  and  voltage  variations  are 
seen  to  be  much  greater  in  Fig.  13  than  in  the  earlier  fig¬ 
ures.  The  error  in  the  algorithm  is  also  larger,  but  the  error 
in  pure  coulomb  counting  is  also  compromised  at  high  rates 
of  charge,  and  it  is  not  unexpected  that  the  combined  SOC 
is  less  than  that  calculated  by  purely  coulomb  counting,  as 
larger  charge  inefficiencies  [27,28]  are  not  included  in  the 
pure  coulomb  counting. 

The  previous  analyses  indicate  that  the  algorithm  is  suf¬ 
ficiently  accurate  of  the  range  of  conditions  investigated. 
Once  the  OCV  and  high-frequency  resistance  R  have  been 
extracted,  it  is  possible  to  estimate  the  charge  and  discharge 
power  capability  of  the  battery  pack,  which  is  important  in 
terms  of  operating  the  vehicle  efficiently  (i.e.,  such  knowl¬ 
edge  is  key  to  determining  optimal  torque  allocation  with 
respect  to  the  heat  engine  and  the  electric  motors).  The  max 
discharge  power  can  be  expressed  as: 


max,  discharge 


=  IV  =  ivmin. 


That  is,  when  the  battery  voltage  obtains  its  lowest  accept¬ 
able  value,  the  max  discharge  power  results.  First  consider 
an  ohmic  battery,  wherein  the  superposition  integral  can  be 
ignored.  For  the  ohmic  battery,  V  =  Voc  —  IR ,  and 


max,  discharge 


=  IVm  in  = 


(Voc  ~  Vmin) 

R 


V™ 


Similarly,  the  max  charge  power  for  the  ohmic  battery  is 
given  by 


max,  charge 


—  TVmax  — 


(Voc  Vmax)  1 

R 


The  ohmic  battery  does  not  address  transient  effects  such 
as  those  correlated  by  the  superposition  integral.  To  improve 
the  estimate  by  including  transient  effects,  we  employ  Eq.  (6) 
to  determine  the  current: 


(Voc  ~  V)t  +  exp(-£  At)[V  -  (Voc  -  IR)]t-At 
R  ±  A  At 


ironic  that  while  the  flatness  in  the  SOC-OCV  curve  gives 
rise  to  great  difficulty  in  determining  the  NiMH  battery  SOC, 
it  is  this  same  lack  of  variation  in  the  OCV  that  allows  NiMH 
batteries  to  provide  nearly  constant  charge  and  discharge 
power  capability  over  a  large  range  of  SOC’s,  which  is  use¬ 
ful  when  it  comes  to  constructing  optimal  vehicle-control 
algorithms. 


5.  Conclusion 

Efficient  vehicle  operation  requires  the  knowledge  of  the 
battery  state;  i.e.,  how  much  energy  is  left  in  the  battery  and 
how  much  power  the  battery  discharge  or  accept.  In  addi¬ 
tion,  on-board  diagnostics  for  the  battery  in  a  hybrid  electric 
vehicle  are  required.  The  algorithm  developed  and  imple¬ 
mented  in  this  work  represents  a  potential  means  to  address 
these  issues.  Calibration  changes  from  values  employed  in 
this  work  are  to  be  expected,  but  the  overall  approach  is 
seen  to  yield  a  substantially  stable  and  robust  algorithm 
over  a  fairly  broad  range  of  application.  More  testing  is 
needed,  however,  to  ascertain  the  errors  in  the  algorithm 
and  conditions  under  which  the  algorithm  might  fail  to  per¬ 
form  acceptably.  Specific  conclusions  from  this  work  are  as 
follows. 

1 .  A  model-reference  adaptive  system  based  on  the  equiv¬ 
alent  circuit  displayed  in  Fig.  4  provides  SOC  accuracy 
to  within  ±10%  for  the  Precept  battery  system  over  the 
Federal  Test  Procedure  schedule.  However,  harsher  drive 
schedules  may  yield  less  accurate  results,  and  the  mea¬ 
surement  error  employed  in  this  work  may  be  of  the  same 
level  of  error  as  that  of  the  algorithm  under  high-power 
excitation. 

2.  NiMH  batteries  exhibit  a  large  hysteresis  voltage,  as  the 
voltage  versus  SOC  curve  under  effectively  equilibrium 
conditions  is  dependent  on  whether  the  battery  is  be¬ 
ing  charged  or  discharged  (cf.  Fig.  3).  This  effect  is 


max,  discharge 


=  /Vmin  = 


(Voc  -  Vmin)r  +  exp(-Z?  A Q[V  -  (Voc  -  /7?)L-a/ 
R  ±  A.  At 


and 

^max,  charge  =  ^Vnax  = 

To  implement  these  equations,  the  respective  powers  are 
calculated  immediately  after  the  algorithm  has  been  em¬ 
ployed  to  finish  the  SOC  determination  at  time  t.  In  this  case, 
quantities  calculated  or  measured  at  time  t  are  then  stored 
in  the  variables  listed  in  the  respective  power  expressions  at 
time  t  —  At.  Then  one  must  state  the  duration  corresponding 
to  the  desired  estimate  for  power.  It  is  sufficiently  accurate 
not  to  update  Voc  unless  At  were  to  exceed  many  seconds. 
One  could  then  use  coulomb  counting  on  the  projected  cur¬ 
rent  to  estimate  a  new  state  of  charge  and  update  Voc.  It  is 


kmax  • 


more  pronounced  at  higher  currents.  We  show  that 
an  adaptive  hysteresis  contribution  is  critical  to  the 
algorithm’s  successful  implementation. 

3.  When  the  algorithm  is  started  arbitrarily  with  an  incor¬ 
rect  initial  state  (e.g.,  incorrect  initial  SOC  and  hystere¬ 
sis  voltage),  about  100  s  are  required  for  convergence.  In 
practice  this  arbitrary  initial  condition  should  not  occur, 
but  the  result  implies  algorithm  stability  and  robustness. 
Such  behavior  should  facilitate  service  procedures  on  the 
battery  pack  as  well. 


(Voc  -  Vmax)r  +  exp (-B  AQ[V  -  (Foc  -  IR)]t _A, 
R  ±  A  At 
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giving  rise  to  the  factor  (It- At  +  4)/2  multiplying  on  to  A. 
Higher-order  approximations  could  be  employed.  Since  the 
expression  is  placed  in  an  adaptive  routine,  and  the  time 
steps  in  practice  are  to  be  taken  as  small  relative  the  char¬ 
acteristic  time  MB,  effectively  identical  results  are  obtained 
with  higher-order  approximations. 
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Appendix  A 


Derivation  of  voltage  recursion  relation .  In  this  section 
we  derive  the  voltage  recursion  relation.  We  begin  with  the 
superposition  integral  solution  of  Eq.  (3): 


V  =  Voc  -  IR  + 


I(Oexp[-B(t  -  0]df. 


This  expression  is  based  on  the  assumption  that  I(t)  =  0 
for  t  <  0.  At  time  zero, 


(V-Voc+IR)t=o  =  0. 

For  the  first  time  step, 

(V  -  Voc  +//?)»=»,  ~  -e~BtlA(It=tl)eBt'  Am, 


or 

(V  -  Voc  +IR)t=tleBt]  ~  -A(/f=fl)eSr>  Am, 

where  the  current  in  the  integrand  has  been  taken  constant 
over  the  time  step.  For  the  second  time  step, 


(V- Voc  +IR)t=t2 

~  e~Btl  [-A(It=tl)zBh-  A *2  -  A(/r=?1)eBf>  Am]  , 
or 


(V  -  Voc  +  IR)l=h 

~  -A{It=t2)At2  +  e~BAt2(V  -  Voc  +  IR)t=h  ■ 

Following  this  procedure,  similar  recursion  relations  re¬ 
sult  for  subsequent  time  steps,  with  the  general  expression 
corresponding  to 

v,  =  (Voc  -  Wt  -  Hit) At  +  e-BAf(V  -  Voc  + 

This  corresponds  to  Eq.  (6).  The  utility  of  the  recursion 
relation  stems  from  the  fact  that  time  integration  of  past 
currents  need  not  be  used  to  determine  the  voltage;  only 
the  present  and  previous  time- step  values  are  required.  It 
should  be  recognized  that  there  are  other  methods  to  approx¬ 
imate  the  integral  that  are  equivalent  from  an  implementa¬ 
tion  point  of  view  when  it  comes  to  constructing  a  recursive 
voltage-current  relation  for  SOC  and  SOH  extraction.  For 
example,  one  can  average  the  current  over  the  time  step, 
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